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ABSTRACT 

This paper is a sequel to our previous work for accretion onto a Schwarzschild 
black hole and the so-called standing accretion shock instability (SASI), in this 
paper we investigate non-axisymmetric perturbations for a Kerr black hole. The 
linear and non-linear phases for the shock evolution are analyzed in detail by both 
2D general relativistic hydrodynamical simulations and linear analysis. Since the 
structure of steady axisymmetric accretion flows with a standing shock wave is 
very sensitive to the inner transonic flow, their properties such as Mach numbers, 
which are important for the stability, depend on the Kerr parameter very much. 
Although the essential features of the instability do not differ from the previous 
results for the Schwarzschild black hole, the frame dragging effects specific to 
the Kerr black hole is also evident. Interestingly, the oscillation periods of the 
fundamental unstable modes are dependent only on the shock radius irrespective 
of the injection parameters. 

Subject headings: black hole physics, hydrodynamics, instabilities, relativity, 
shock waves 



1. Introduction 

It is widely believed that the transonic rotating accretion flow onto a black hole is 
one of the candidate sources for many high energy astrophysical phenomena. The time 
variability of these objects may be interpreted as the unsteadiness of accretion flows, which 
in general suffer from various types of instabilities, such as thermal, viscous, dynamical, self- 



gravity, induced ones and so on. A variety of physical mechanisms (IVarnie fc Taggerl (12002 
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20061 ); iDong fc David I (120081 ). and references theirin) have been considered to explain the 
observed intensity variations. The shock wave in the accretion flows is one of the promising 
mechani sms and the so-called "shock oscillation model" has been investigated by ma m 



mecnam sms and tne so-called snock oscillation model nas been investigated by ma ny 
authors flchakrabartj[l989l: IMolteni et alilrgyj^akrabarti et.al.l booi bas et al.l l2003al lbl: 



Aoki et al. 



2004 



Molteni et al 



20061 ; lOkuda et al.l 120071 ). The stability of accretion flows 



with a standing shock wave in them is one of the key issues. 

The linear stability of the standing shock wave against axisymmetric perturbations in 
the axisymmetric accretion flow is well-known and simply stated a s follo ws: if the post- 
shock matter is accelerated, the shock is unstable (INakayamal 1 19941 . Il995l ). However, the 
accretion shock wave that is stab le against axisymme t ric perturbations can be unstable for 



non-axisymmetric perturba tions (IMolteni et al.l 119991 ; iGu fc Foglizzol 120031 ; iGu fc Lul 12005 
Nagakura fc Yamadal 120081 ) . This non-radial shock instability may be quite common. In 
fact, the similar instability is found in accretion flows onto a nascent proto neutron star in 
the supernova explosion and might be an important element in the mechanism for it. 



Our previous paper (INagakura fc Yamadal (120081 ) hereafter NYvl) investigated general 
relativistically the non-axisymmetric shock instability around a non-rotating black hole, 
suggesting also possible implications for fluctuations in GRB jets. However, there is no doubt 
realistic stellar mass black h oles in general are rotating as suggested by the observations of 
massive stars (IFukudal Il982l ). Since it is believed that rapidly rotating massive stars are 
the progenitor of long GRBs, the investigation of the standing accretion shock instability or 
SASI should be extended to the Kerr black hole. Moreover, in order to take into account 
the frame dragging due to the black hole's spin, the fully general relativistic treatment is 
strongly required. We observed in the previous paper that the general relativistic effects are 
important to determine the structure of the inner transonic flow. This in turn means that 
the injection parameters such as Bernoulli constant and specific angular momentum are also 
affected by the general relativistic effects. In fact, the injection parameters that allow the 
existence of the standing shock wave are sensitive to the Kerr parameters. 

The plan of this paper is as follows. In section 2, we describe the steady axisymmetric 
accretion flows with a standing shock wave in them. Then we present the formulation of 
linear analysis in section 3. The numerical methods for dynamical simulations and models 
are described in section 4. The results of the numerical simulations obtained by both linear 
analysis and dynamical simulations are shown in section 5. We conclude the paper in section 
6 with summary and discussion. 
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2. Axisymmetric steady accretion flows with a shock wave in Kerr spacetimes 



The mathematical formalisms of axisymmetric steady accretion flows onto Kerr black 
hole with a standing shock wave are essentially the same as those used in NYvl. Only the 
equatorial plane is considered and the Kerr space metric is fixed in time. The conditions 
for the existence of multi ple sonic points and sta nding shock waves in Kerr spacetimes 



were already discussed by (jLulll986 



Lu et al.l 119971 ). In the following we briefly review the 



essential points for constructing axisymmetric steady flows , whic h would he l p the reader to 
understand the property of the solutions. We refer to|Lu| (119861 ); iLu et al.l (119971 ) for more 
details. 

We first determine the range of the injection parameters such as Bernoulli constant and 
specific angular momentum, these give multiple sonic points. The existence of multiple sonic 
point is a necessary condition for the existence of standing shock wave, since not only the 
pre-shock flows but also post-shock flows are transonic. Then, two transonic accretion flows 
are connected by a standing shock wave, where the Rankine-Hugoniot relations are satisfied. 
Since the inner transonic flows are strongly affected general relativistic effects, the structure 
of the entire accretion flows also varies with Kerr parameters. 

The basic equations are the relativistic continuity equation and equation of energy- 
momentum conservations. 



0. 
0. 



(1) 

(2) 



Po,u fJ, ,T' lu are the rest mass density, 4-velocity and energy momentum tensor, respectively. 
Throughout this paper, we assume a T-law EOS, p — (T — 1) p e, where p and e are the 
pressure, and specific internal energy. Except for the discussion of dynamical simulations, 
we adopt the geometrical units G = c = 1, where G and c are the gravitational constant 
and the speed of light. Latin and Greek indices denote spatial and spacetime components, 
respectively. Since we consider the accretion flows only in the equatorial plane, we ignore 
the ^-component of velocity, u , and all 9 derivatives. The basic equations can be reduced 
to ordinally differential equations with respect to the radial coordinate: 



d r (r 2 p u r ) 
d r p + pou r d r (hu r ) 

d r (hu t ) 

d r (hlic/,) 
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Po h { (d r g rr ) {u'f + {d r g H ) (O 2 + (d r g tt ) (uf + 



2 (d r g t t) u l u* 
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where we use the Boyer-Lindquist coordinates. Here, g^, h and p are the metric, spe- 
cific enthalpy and pressure, respectively. The only difference from the equations for the 
Schwarzschild black hole is the last term in eq. (jlj), {2 (d r g t( p) which represents the 

frame dragging effect. It should be noted, however, that the frame dragging effect is included 
not only in this term but also in the 4-velocity itself. The tran s onic fl ow can be obtained by 



imposing the regularity conditions at the sonic point (see |Lu| (119861 ) ). where the equations 
become apparently singular. Figure [I] shows the locations of the sonic points as a function 
of the specific angular momentum and Kerr parameter. The adiabatic index and Bernoulli 
constant are fixed as 4/3 and 1.004, respectively in this figure. As is clearly demonstrated, 
the curve is symmetric with respect to A = 0, which A denotes the specific angular mo- 
mentum, for the Schwarzschild spacetimes, where as it is not true of the Kerr spacetimes. 
This is due to the frame dragging of Kerr black holes, which acts effectively as an additional 
angular momentum that forces a corotation of matter with the black hole, and thus the 
symmetry is broken. This figure also shows that there are indeed multiple sonic points for 
adequate injection parameters. It is well known that the physically-meaningful sonic points 
are the inner and outermost ones, which are so-called saddle-type singularities, while the 
unphysical middle sonic point is a central-type singularity. Hence, the steady accretion flows 
with a standing shock wave should pass through the inner- and outermost sonic points. The 
Rankine-Hugoniot relation determines the location of standing shock wave from these two 
transonic flows. These equations are completely the same as those used in the previous paper 
(see Section 2.2 in NYvl). Just like sonic points, the shock positions also vary with the Kerr 
parameter very much. 

As in the Schwarzschild black holes, there are generally two possible shock locations 
under adequate injection parameters. These are referred to as the inner and outer shock 
points. We consider only the outer shock in this paper, since it is well known that the inner 
shock is alr eady unstable aga inst axisymmetric perturbations, which is also true of the Kerr 



spacetimes (jNakayamal Il995l ) . We construct axisymmetric steady accretion flows imbedding 
on outer shock wave for different injection and Kerr parameters, which are summarized in 
Table [H The adiabatic index is fixed to T = 4/3 throughout the paper. It should be noted 
that some models have negative specific angular momentum, i.e, matter rotates in the oppo- 
site direction of the black hole spin, although it is rather unlikely in realistic astrophysical 
systems. It is, however, useful to elucidate the Kerr parameter dependence and the following 
fact, which will be demonstrated shortly, also gives another merit: as the Kerr parameter 
gets larger, the range of the injection parameters that allows the existence of a standing 
shock wave become smaller and disappears eventually in the corotation case, whereas it is 
not true for the counter rotation case. This allows us to investigate the properties of the 
instability under wider injection parameters. 
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In Figure [21 we show the ranges of injection parameters, for which the standing shock 
wave can exist. As seen in the left panel, where the corotation case is presented, as the 
Kerr parameter becomes larger, the allowed regime of the specific angular momentum shifts 
to smaller values. This is due to the frame- dragging by the Kerr black hole, which acts 
effectively to increase the angular momentum of matter. Note that the specific angular 
momentum is conserved along the stream line and frame- dragging effect is not restricted to 
the vicinity of black hole. The boundary of the allowed regions normally consists of two 
areas. For a given Bernoulli constant, the maximum angular momentum is determined by 
the weak shock limit, i.e, Mach number of the shock M s h ~ 1. On the other hand, the 
minimum value is determined by the location of the outer shock wave. It is widely known 
that the outer shock wave is located at a larger radius than the middle sonic point, while the 
inner shock wave is located at a smaller radius. If the specific angular momentum becomes 
smaller, the locations of each shock wave approach the middle sonic point, and eventually 
coincide at this point. This determines the minimum specific angular momentum for the 
existence of a standing shock wave. As mentioned in NYvl, however, for small Bernoulli 
constants, there is a maximum specific angular momentum above which no multiple sonic 
points exist. This is true only the fully general relativistic treatment. Hence, there is a 
possibility that the maximum specific angular momentum is determined by this limit. As 
a matter of fact, there is indeed the case of a/M* = 0.9, in which the allowed regime is 
bounded by a vertical like for small Bernoulli constants. 

The right Figure [2] shows the region of injection parameters, for which the steady shock 
wave exists in the conter-rotation case, where matter rotates in the opposite direction of the 
Kerr black hole spin. As expected, the allowed regions are quite different from those for the 
corotation case. As the Kerr parameter becomes large, the specific angular momentum gets 
more negative and the Bernoulli constant becomes smaller. Because the frame dragging ef- 
fects decelerate the matter rotation, the specific angular momentum has to be more negative. 
As mentioned above, the specific angular momentum is conserved along the stream line, it 
has too large an absolute value (and hence too great centrifugal force) at distant places from 
the black hole, where the frame dragging is negligible. The excess should be compensated 
for by reducing the Bernoulli constant and as a result, pressure. 



3. Linearized equations in Kerr spacetimes 

We perform a linear analysis for the non-axisymmetric shock instability around a Kerr 
black hole. The obtained eigen-values help us to analyze not only the shock stability, but also 
validate the accuracy of dynamical simulations. By combining dynamical simulations with 
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linear analysis, both the linear and non-linear phases can be investigated in detail. In this 
section, we present a formulation for the linear analysis including the boundary conditions. 

The basic equations are the linearized relativistic continuity and energy-momentum con- 
servation equations. The perturbed quantities are assumed to be proportional to e - lu)t + im <t> . 
Since we consider the equatorial plane only and neglect all the ^-derivatives and the 9- 
component of velocity, the linearized equations are expressed as follows: 

drf = '—^{poa^ci + po^iuu^-mu^)}, (7) 

d r q = pmU r ( o )h{0)Um ("PW+Po(o)u m h {0) u m a q ) , (8) 
u t(o) 

u t(0) 

d r S {1) = i—^aSp.), (10) 
where S is the entropy per baryon and the following notations are used: 

/ - ^ + (») 

Po(o) u \ * 

» = ^>+^, (12) 

tyo) «t(0) 

V(i) = to ( fyi) 1X^(0) + ^(0)^(1) ) + m ( h(i)U t (p) + h(Q)U ti x) ) , (13) 

/ and q are the values which relate to the mass flux and Bernoulli constant, respectively. V\ 
is a useful value which we can integrate analytically. The a denotes the complex frequency 
of the perturbations measured in the rotating frame. Note that these equations for the 
Kerr spacetimes look completely the same as those used in NYvl for the Schwarzschild 
spacetimes. As mentioned in the previous section, however, the frame-dragging effects are 
included implicitly in the components of 4- velocity. As in the Schwarzschild black hole case, 
V(i) and Sn) can be integrated analytically (see Section 3 in NYvl), and we only need to 
integrate numerically Eqs.(j7J) and (jSJ). 

The boundary conditions are also the same as those used in NYvl. We impose the outer 
boundary condition on the linearized Rankine-Hugoniot relation at the shock surface, assum- 
ing that the pre-shock region remains unperturbed. The explicit forms of these equations 
are given by, 



(p (o)^ (1) ) + +(po(i)^ {0) ) + + ^ = 0, (15) 
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(16) 



1 VPo(o) tyo) u m u r ^>J) + 

{(MurWf* + ^ + ^ + ^) +«.,} + C = 0, (17) 

>> VP0(0) fyo) M r(0) M r( °7 J + 

{ (Mtv0 <°> ( m + M + m + ^) } + D = o, us) 

>> VPo(o) fyo) «#o) M r(0) / J + 

with the following definitions: 

A EE ^{( Po(0) ^ (19) 

£ ee iwrj | (pohutu 1 + p)| 0) - (pohutu* + p) W j 

-imr) { (p o /w ) f - (po/w )^ } , (20) 

+i<^?7 1 (p ^*-u r )| 0) - (po^*-u r )^ 0:i | - im?] | (p /iM r -u^)| 0) - (p /i-u r -u^) (21) 
£ ee ia;?7 { (po/wt^,?/) j 0) ~~ (po^V"*)^} 

{ (PO^^ + P)| 0) - (p + P)i 0) } • (22) 

where 77 denotes the displacement of the shock radius and defined as: 

Rsh = R s h(o) + V exp ( -iotf + imcj) ) , (23) 

where R s h stands for the shock radius. 

Finally, the regularity condition is imposed at the inner sonic point as the inner boundary 
condition. The explicit form of this condition is; 

G ( i) - F ( i)U r( o), r = 0, (24) 

where the following notations are employed 

F (0) ee p (o)V)/ r {M r( V(o) - (& S (0)) 2 (1 +« r{0) M r(0) )} (25) 
F (1) ee p 0(1) h {0) (u r ^) 2 + Pm h (l) (u^f + 2p m h {0) u r ^u r ^ 

-rg rr [ P{1) (1 + M r( V(o)) + 2p (0) u r(1) u r (o)}, (26) 

G(i) EE T f 5( rr ir + - J I (p(l)ttr(O) +P(0)«r(l)) (l + Wr(0)« r(0) ) + 2p (0 ) (u r ( )) 2 M r(1) } 
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+ ^{P0(l)^(0)w r(0) + P0(0)fyl)M r(0) + P0(0)fy0)W r(1) 
{ 9 rr,r (u^f + 9 H , (u^f + 9tt , («*<°>) 2 + 2 9t4> ,U^U^} 

+ iTp (0 ) (1 + u r(0) u r(0) ) (-cou t(1) + mu^ 1] ) 

-iu m P(l) a. (27) 

The unperturbed sound velocity in the comoving frame and the adiabatic index are de- 
noted by 6 S ( ) and T, respectively. Note again that the only explicit modifications from the 
Schwarzschild spacetime case are the terms with g t( p (see section 3 in NYvl), although the 
frame dragging effects are included implicitly in the 4-velocity. 



4. Numerical Methods and Models 

We perform two dimensional general relativistic hydro dynamical simulations of non- 
axisymmetric shock instability in the equatorial plane around a Kerr black hole. The nu- 
merical methods in the present paper are essentially the same as those used in NYvl. 

Our numerical code is based on the so-called central scheme, which guarantees a good 
accuracy even if flows include shock waves. We refer the reader to the Appendix of NYvl 
for more details on this code. The Kerr-Schield coordinates, which have no coordinate 
singularity and allow us to put the inner boundary inside the event horizon, are used. The 
computational domain of 0.6r eve <r< 200M* is covered by the 600(r) x 60(0) grid points, 
where M* and r eve stand for the black hole mass and the radius of the event horizon in the 
equatorial plane, respectively. The radial grid width is non-uniform, increasing by 0.34% 
per zone from the smallest grid (Ar = 0.1M*) at the inner boundary. The computational 
times are ~ 6 x 10 M*, which corresponds to ~ 500ms for M* = 3M Q . In the following, we 
discuss the numerical results in cgs units for the case M* = 3M , since we have in mind the 
application to massive stellar collapses as in NYvl. Note, however, that the formulations 
are dimensionless, and the scaling is obvious. 

The initial perturbation mode has m = 1, where m indicates the azimuthal mode number 
in e tm ^, and the amplitude is fixed to 1% (see NYvl for the dependence on the mode number 
and the initial amplitude). The radial distributions of initial perturbations are obtained by 
the linear analysis. We choose the most unstable m = 1 mode, which allows us to compare 
clearly the linear growth phase between the linear analysis and dynamical simulations. 
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5. 



Numerical Results 



In this section, we first describe the results of the linear analysis, which are obtained 
by solving Eqs. fl7|)- f|T0|) with the appropriate boundary conditions at the shock surface 
and sonic point (see Section 3). Then, we present the time evolutions of shock instability 
obtained by GRHD simulations in both linear and non-linear phases compared them with 
the linear analysis. In so doing, we employ the mode decomposition on the shock surface by 
the Fourier transform: 



Jo 

where R s h {<p, t) and a m (t) are the radius of the shock wave as a function of </> and t and the 
amplitude of mode m as a function of t, respectively. 



The real and imaginary parts of unstable eigen-frequencies for the standard Model 
Mk05E4L301p are given in Figure [31 which also shows the another parts for Model Mk00E4L343p 
for comparison. In these models, the Bernoulli constant is identical and the Kerr parameter 
and specific angular momentum are chosen so that the shock radius would become the same. 
As seen in Figure |3j there are the unstable eigen-frequencies for Model Mk05E4L301p, and 
we also find that the higher frequency is stable than the lower frequency among these models. 
It should be noted, however, that this property is not always satisfied (see Fig. 5 in NYvl). 
The magnitude of the imaginary part of omega is slightly larger than the Mk00E4L343p 
case, because the shock wave for Mk05E4L301p is stronger than Mk00E4L343 case. 

Figure [3] also tells us the dependence on the Kerr parameter: the flow tends to be 
more unstable as the Kerr parameter gets larger. Note that the Bernoulli constant is not 
changed and the specific angular momentum is chosen in such a way that the shock location 
is unchanged. This trend is also seen in other models. In Figure HI we show that the real and 
imaginary parts of unstable eigen-frequencies for models Mk00E4L350p to Mk0999E4L406i 
(See Tabled]), have the same shock position (R sh = 35M*). As can be clearly seen in this 
figure, the imaginary part of the eigen-frequency becomes larger in the corotation case, while 
it becomes smaller in the inverse rotation case (with increasing Kerr parameters). The reason 
is as follows. 

From Table [U we find that the growth rates and the Mach numbers of the pre-shock 
flows are correlated with each other. Stronger shocks tend to have a larger imaginary part 
of eigen-frequency, i.e., they are more unstable. This fact is the key to the above trend. The 




(28) 



5.1. Linear Analysis 
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inner transonic flow is strongly affected by the frame dragging of the Kerr black hole and it 
works effectively as an addition to the specific angular momentum of matter. As a result, the 
specific angular momentum becomes smaller with increasing Kerr parameters if the shock 
location is fixed. Since the outer transonic flow is hardly affected by the frame dragging 
effect, this leads to the reduction of the centrifugal force there. As a result, the radial 
velocity of outer transonic flow becomes faster (slower) in the corotation (counter-rotation) 
case, and the Mach number of the pre-shock flow gets larger (smaller) with increasing Kerr 
parameters, which then leads to the trend that the models with larger Kerr parameters tend 
to be more unstable (stable) in the corotation (counter-rotation) case. 

It is also interesting to present the dependence of shock stability on the Kerr parameter 
under the same injection parameters. In Figure \5[ we shows the eigen-frequencies for models 
Mk00E4L345p, Mk005E4L345p and Mk01E4L345p, all corotation cases. The initial steady 
flows for these models are constructed with the same injection parameters but different Kerr 
parameters. As found from Table [TJ the shock radius becomes larger with increasing Kerr 
parameters in this case. Recalling our previous results for the Schwarzschild black hole, 
that the shock radius becomes larger as the specific angular momentum is increased with 
the same Bernoulli constant, we can interpret the present results as follows: as the Kerr 
parameter increases, the specific angular momentum becomes effectively larger, and the 
shock is pushed outwards. Then the Mach number of pre-shock flows becomes smaller, since 
the radial velocity is a decreasing function of radius. As a result, the flow becomes more 
stable, or Ui gets smaller. 

In Figures H] and [5j we find that the unstable eigen-frequencies for higher harmonics 
oscillate for some models. We think that this is not a numerical ar tifact, and the oscillatory 



behavior is found in preceding papers (NYvl. lGu fc Foglizzd (120031 )). At present, the reason 



for these oscillations is unclear, and further investigations are necessary. 

Finally, we discuss a possible strong correlation between u r for the fundamental unstable 
modes and the shock radius. In Figure El we show the shock radii and u r for the fundamental 
unstable modes for all the various models. We also include the previous results (NYvl) in 
this figure to increase the number of samples. The injection parameters for the previous 
results are listed in Table 1 in NYvl. As is evident, the real part of fundamental unstable 
modes decreases monotonically with increasing shock radius even though all models does not 
have the same injection parameters, adiabatic indices or Kerr parameters. This suggests that 
the shock radius is the unique factor to decide the oscillation frequency. It also means that 
the rotational pattern velocity of the deformed shock surface, which corresponds to u r , is a 
function of the shock radii alone. It is interesting that the specific angular momentum has no 
direct effect on it. The relation will also be useful to infer u r from the shock radius without 
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detailed calculations. At present, we do not understand the reason for this correlation. This 
may be a clue to the understanding of the mechanism SASI, and definitely needs further 
investigations. 

5.2. Numerical Simulations 

Although we investigate several models with different Kerr parameters, the qualita- 
tive evolutionary path is generally similar to that found in the previous simulations NYvl. 
We first summarize the basic features in the temporal evolution of the standard Model 
Mk05E4L301p, which has the same Bernoulli constant and shock location but a Kerr pa- 
rameter and specific angular momentum different from Model Mk00E4L343p, which is the 
standard model Ml in the previous paper (NYvl). Figure [7] shows the snapshots of the 
velocity for this model. In the linear phase, a purely single m = 1 mode grows exponentially 
and the deformed shock wave starts to form a spiral arm, which rotates in the same direction 
as the unperturbed flow. Over times, the spiral arm develops and non-linear mode couplings 
can not be negligible. When the non-linear regime is reached, the dynamics becomes highly 
complex. Several shocks start to collide with each other and the original shock begins to 
oscillate radially. While large and small radial oscillations are repeated like a limit cycle, we 
find that quasi-steady of asymmetric configurations m = 1 or 2 are realized eventually for 
most of the models. For some models that have too large or too small radial oscillations, 
such as models Mk01ElL350p and Mk02E8L350i, we can not identify the dominant mode 
in the non-linear phase. 

In Table H] we give the saturation amplitudes of m = modes, which vary from model 
to model. We recognize a strong correlation between the Mach number and the saturation 
amplitude of axisymmetric m = mode: stronger shocks tend to give larger saturation 
amplitudes. As a matter of fact, the shock wave for Model Mk01ElL350p, which was the 
largest Mach number, leaves the computational domain soon after the beginning of numerical 
simulations. On the contrary, the shock wave for Model Mk02E8L350i with the smallest 
Mach number does not move so much from the initial position throughout the simulation 
(see Figured]). 

The oscillation timescale and growth rates obtained from the dynamical simulations 
are in good agreement with those found in the linear analysis (see the left panels of Fig- 
ures M and [10]) . Thus, we think that the dynamical simulations are accurate enough. As a 
matter of fact, we have also done the simulations with doubled spatial resolutions both in 
the r and <fr directions. The computational domain is covered by 1200(r) x 120(0) grid points 
with the smallest grid of Ar = 0.05M* at the inner boundary. The difference is roughly 10% 
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at most (see Figure [XT]) . 

It is also found from these figures that the results of dynamical simulations start to 
deviate from those expected by the linear analysis after ~ 6ms for Model Mk05E4L301p and 
~ 30ms for Model Mk0999E4L406i, which indicates the beginning of the non-linear phase. 
In fact, in the right panels of these figures, which show the evolution of the axisymmetric 
m = mode, we find that it starts to grow from the time at which the non-linear phase 
begins. It is important to point out that the duration of the linear phase is much longer in 
Model Mk0999E4L406i than in Model Mk05E4L301p, since the former is more stable than 
the latter. 

Finally, we compare the obtained growth rates with the acoustic-acoustic and the 
advective- acoustic cycles between the shock surface and the inner reflection point. This 
is done to shed light on the instability mechanism. In this analysis, we assume that the 
inner reflection point is either the corotation radius or the inner sonic point. The procedure 
in the comparison is completely the same as that used in NYvl (see Section 5.5 in NYvl). 
The ratios of the growth rate to the period of each cycle are presented in Figure [TSJ It 
is observed that the instability is more likely to be of the Papaloizou-Pringle type. Note, 
however, that in general the wavelengths of acoustic perturbations are longer than the scale 
heights (see Table [2]) and we can not clearly identify the reflection point, just as in NYvl. 



6. Summary and Discussions 

This paper is an extension of our previous work NYvl. The non-axisymmetric standing 
accretion shock instability around a Kerr black hole has been investigated by performing 
both linear analysis and dynamical simulations with fully general relativistic hydrodynamical 
treatments. This allows us not only to make clear the dependence of SASI on the Kerr 
parameter, but also to elucidate the effects specific to the Kerr black hole, such as the frame 
dragging. 

We have explored the range of the injection parameters that allows the existence of a 
standing shock wave in axisymmetric steady accretion flows and made clear its dependence on 
the Kerr parameter. In particular, we have demonstrated for large Kerr parameters and small 
Bernoulli constant in the corotation case, that the maximum specific angular momentum for 
the existence of standing shock wave is determined by the limit for the multiple sonic points. 
This is qualitatively different from the limit obtained from the pseudo-Newtonian analysis 



f lGu fc Foglizzcl 12003 ) 



According to the linear analysis, the shock wave becomes unstable against non-axisymmetric 
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perturbations, and the frame-dragging effects are also evident. The growth rate of the in- 
stability varies widely with the Kerr parameter. Since the inner transonic accretion flow is 
strongly affected by the frame-dragging from Kerr black hole, the specific angular momen- 
tum of matter near the black hole gets effectively increases. If the shock radius is unchanged, 
the intrinsic specific angular momentum of matter becomes smaller as the Kerr parameter 
gets larger. As a result, the Mach number of the pre-shock flow becomes larger (smaller) 
with the increasing Kerr parameter for the corotation (counter-rotation) flow. We have also 
investigated the dependence of the instability on the Kerr parameter for the same injection 
parameters instead of the same shock radius in the corotation case. The location of standing 
shock wave is shifted outward as the Kerr parameter increases, since the inner transonic flow 
gains effectively the specific angular momentum from the frame- dragging by the black hole, 
and the centrifugal force is enhanced. We have found that the oscillation frequency of the 
non-axisymmetric shock wave depends only on the shock radius, irrespective of the injection 
parameters as well as the Kerr parameters. At the same time, we find a correlation between 
uj r for the fundamental unstable mode and the shock radius. The oscillation frequency de- 
creases monotonically with increasing shock radius. Thus, we conclude that the shock radius 
is the unique factor to decide the oscillation frequency. 

In the dynamical simulations, we have observed that the essential evolution of the 
perturbed system does not change so much from that for the Schwarzschild spacetimes. We 
have seen that the quasi-steady configuration of m = 1 or 2 nature are formed eventually 
for most of the models, which were also observed in NYvl. The maximum saturation level 
of average shock radius has been found to be correlated with the Mach number of the pre- 
shock flow, which is correlated in turn with the eigen-frequency obtained from linear analysis. 
Hence, the dependence of the saturation level in the non-linear phase on the Kerr parameter 
can be understood by the result of linear analysis. 

It is important to point out that although the ranges of the injection parameters that 
allows the existence of st anding shock waye is not very wide, they are obtained indeed in 



recent progenitor models (IHeger et al.ll2005l ). In fact, if we assume that the black hole mass 
is M* = 3M and the temperature of matter is less than 1 MeV, the range of specific angular 
momentum for the existence of shock wave is about X cgs ~ 3 — 5 x 10 16 cm 2 /s. Since the 
standing shock wave is very robust, we believe that a shock wave is formed one way or another 
provided the injection parameters are appropriate. If such an accretion flow produces GRB 
somehow, the SASI will be a natural source of fluctuations in the GRB jets. 

In order to apply the present results to realistic astrophysical phenomena, such as GRB 
and QPO, there are many improvements remaining to be made. The radiation, viscosity and 
magnetic field should be investigated in addition to making the model more self-consistent, 
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all of which are currently under way. Moreover, needless to say, the structure in the meridian 
section should be taken into account, which has been a major challenge for many years, since 
the accretion flow on to the black hole is inevitably transonic and it is difficult to impose a 
regularity condition at the sonic surface. In t his respect, i t is encouraging that the problem 



was solved by the linear perturbation method (lBeskinll2002l ). It should be noted, however, the 



unperturbed flow is assumed to be spherically symmetric and the specific angular momentum 
and the frame-dragging by the Kerr black hole are small in the paper. We will need to obtain 
steady solutions for wider injection and Kerr parameters, which is our future work. 

This work was partially supported by the Grant-in- Aid for the 21st century COE pro- 
gram "Holistic Research and Education Center for Physics of Self-organizing Systems" of 
Waseda University and for Scientific Research of the Ministry of Education, Science, Sports 
and Culture of Japan (17540267, 14079202). 
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Fig. 1. — The locations of sonic points as a function of the specific angular momentum and 
Kerr parameter. The adiabatic index and Bernoulli constant are fixed as 3/4 and 1.004, 
respectively. 
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Fig. 2. — The ranges of the injection parameters that allow the existence of a standing 
accretion shock wave around Kerr black holes. The red, green, blue and pink lines show the 
boundaries of the allowed regions for a/M* = 0, 0.3, 0.6 and 0.9, respectively. The left panel 
shows the corotation case and the right panel gives the counter-rotation case. 
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Fig. 3. — The eigen-frequencies of unstable m = 1 mode for models Mk05E4L301p (stan- 
dard model) and Mk00E4L343p. These models have the same shock location (r = 16M*) 
and Bernoulli constant (E = 1.004), but a different Kerr parameter and specific angular 
momentum. 
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Fig. 4. — The same as Figure [3j but for different models. The shock location of these models 
is fixed to r = 35M*. Note that the negative values of u r correspond to the confer- rotating 
models. 
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Fig. 5. — The eigen-frequencies of unstable models for the models with the same injection 
parameters, but different Kerr parameters. 
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Fig. 6. — The shock radii and the real parts of the eigen-frequencies for the fundamental 
unstable modes for various models. In addition to the present results (green dots), the 
previous results (NYvl) are included as red dots. 
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Fig. 7. — The time evolution of the velocity for Model Mk05E4L301p. The color contour 
shows the magnitude of radial velocity. The arrows represent the velocities at their positions. 
The central region in blue corresponds to the black hole. 
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Fig. 8. — The evolutions of the amplitude of m = mode for Mk02E8L350i which we evolve 
long time (1000ms) simulations. The same evolutions for standard Model Mk05E4L301p are 
also shown for comparison. 
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Fig. 9. — The time evolution for Model Mk05E4L301p. Left: the comparison of the time 
evolutions of the amplitudes of the m = 1 mode between the linear analysis (red lines) and 
dynamical systems (green crosses). Right: the evolutions of the amplitude of the m = 
mode. 




Fig. 10.— The same as Figure M but for Model Mk0999E4L406i. 
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Fig. 11. — The time evolutions of m = 1 (left) and m = (right) modes for the different 
resolution for the standard model Mk05E4L301p. The red line is the result for the standard 
resolution in this paper, the green line is the result for the doubled spatial resolution. 
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Fig. 12. — The ratio of the growth rate to the frequencies of advective-acoustic (+) and 
acoustic-acoustic cycle (x) for all the models. In the left (right) panel we assume that the 
corotation (inner sonic) point is the inner reflection point. See also Figure 10 in NYvl. 



Table 1. Model Parameters 





Kerr Parameter 


Bernoulli 


Specific Angular 


Inner Sonic 


Shock Point 


Mo.cn Number 


Model 


a[M*] 


Constant E 


Momentum A [M*] 


Point ^ c r . \ TV/* 1 


' s/i L * J 




Mk00E4L343p 





1.004 


3.43 


5.3 


16.1 


2.444 


Mk05E4L301p 


0.5 


1.004 


3.01 


3.6 


16.0 


2.682 


Mk00E4T,350r> 


o 


1.004 


3.50 


5.0 


34.8 


2.078 


Mk03E4L327p 


0.3 


1.004 


3.27 


4.1 


35.0 


2.137 


Mk06E4L297p 


0.6 


1.004 


2.97 


3.1 


35.0 


2.208 


Mk03E4L370i 


0.3 


1.004 


-3.70 


5.9 


35.0 


2.019 


Mk06E4L386i 


0.6 


1.004 


-3.86 


6.8 


35.0 


1.967 


Mk0999E4L406i 


0.999 


1.004 


-4.06 


7.9 


35.0 


1.904 


Mk01ElL350p 


0.1 


1.001 


3.50 


4.5 


35.0 


3.426 


Mk01E6L350i 


0.1 


1.006 


-3.50 


5.6 


35.0 


1.632 


Mk02E8L350i 


0.2 


1.008 


-3.50 


6.1 


35.0 


1.395 


Mk00E4L345p 





1.004 


3.45 


5.2 


20.8 


2.355 


Mk005E4L345p 


0.05 


1.004 


3.45 


4.9 


29.9 


2.180 


Mk01E4L345p 


0.1 


1.004 


3.45 


4.6 


44.8 


1.931 



Note. — Parameters for the initial axisymmetric steady flow with a standing shock wave. M* is the mass of the central black 
hole. 



Table 2. Properties of Instability 





maximum amplitude 


Oscillation Period 


Growth time 


for 


Wavelength of Acoustic Perturbations 


Model 


of m = mode 


tosci 


tgrow/2lT 


fundamental node 




Mk00E4L343p 


3.9 


3.7ms 


3.2ms 


1677 


142.6km (31. 7M*) 


Mk05E4L301p 


9.7 


3.7ms 


2.5ms 


1681 


145.8km (32.9M*) 


Mk00-h4L35Up 


3.2 


5.9ms 


8.2ms 


525 


171.9km (38. 2M*) 


Mk03E4L327p 


4.5 


6.1ms 


7.5ms 


517 


174.6km (39.4M,) 


Mk06E4L297p 


5.5 


6.2ms 


6.8ms 


503 


181.3km (40.9M.) 


Mk03E4L370i 


3.3 


5.9ms 


9.3ms 


531 


258.0km (58.2M*) 


Mk06E4L386i 


2.2 


5.9ms 


10.7ms 


535 


244.8km (55.3M.) 


Mk0999E4L406i 


1.8 


11.7ms 


12.4ms 


537 


458.3km (103.4M*) 


Mk01ElL350p 


5.5 


12.5ms 


6.0ms 


504 


343.7km (77.6M*) 


Mk01E6L350i 


1.3 


6.0ms 


14.1ms 


526 


279.1km (63.0M.) 


Mk02E8L350i 


1.02 


4.3ms 


35.7ms 


511 


202.2km (45.6M.) 


Mk00E4L345p 


4.5 


5.4ms 


4.2ms 


1156 


186.5km (42.1M.) 


Mk005E4L345p 


5.6 


4.77ms 


6.7ms 


667 


143.9km (32.5M.) 


Mk01E4L345p 


3.8 


8.67ms 


11.4ms 


357 


225.2km (50.8M.) 



Note. — t osc i, t grow /2iT and \ w represent the oscillation period, growth time and wavelength of acoustic perturbations, respectively, which 
are obtained by linear analysis. We refer the reader to the NYvl for the mathematical definitions of these values. u r (f) i s t ne real part of 
unstable fundamental eigen-frequency for each model. 



